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I. INTRODUCTION 

There have been remarkable experimental advance- 
ments in the area of trapped one-dimensional cold atom 
systems pj, [2, B 0] ■ The one-dimensional quantum gas 
can be realized experimentally by tightly confining the 
atomic cloud in the radial directions and weakly confin- 
ing it along the axial direction la, J™. Since the radial 
degrees of freedom are frozen 0, lf|, the quantum gas 
can be effectively described by a one-dimensional model 
along the axial direction. Using Feshbach resonance, the 
energy of a bound state of two colliding atoms in a mag- 
netic field was tuned to vary the scattering length from 
—oo to oo, which allows experimental access to strongly 
interacting regime and studying the crossover from a BCS 
to Bose-Einstein Condensation. [3, [I(J, [HI, LL2J • The ex- 
perimental achievements have opened up many exciting 
possibilities for theoretical investigation of quantum ef- 
fects in low-dimensional many-body systems. 

Exactly solvable models provide us an important in- 
sight in studying the quantum many-body physics be- 
yond various approximation schemes. The exact solu- 
tions can supply some believable results thus serve as a 
very good starting point to understand the new phenom- 
ena and new quantum states in trapped cold atomic sys- 
tems. Recently, the study of one-dimensional bosons with 
contact interactions in an infinite potential well shows 
a revival of interest [13, LLJ]- We note that the point- 
interacting two-component Fermi gas in a hard wall trap 
is also interesting. Many new physics appears in the 
fermionic system and most important is that this system 
is integrable [U, E ES El El- Meanwhile, this sys- 
tems can be realized in experiment by the laser cooling 
technique. 

In this paper, we investigate the one-dimensional ultra- 
cold two-component Fermi gas in an infinite potential 
well. Most of the previous studies focus on the en- 
ergy spectrum and its related quantities. Here, we put 
our attentions on the wave function of the system. We 
study the single-particle reduced density matrix, natural 
orbitals, the momentum density distributions and two- 



body correlations based on the exact wave function of 
the system. These quantities are important for a bet- 
ter understanding of the ground state properties of the 
fermion system and can be detected in experiments di- 
rectly. 

The paper is organized as follows. In Sec. II, the model 
Hamiltonian and its Bcthe-ansatz solution are introduced 
in details. In Sec. Ill, we study the one-body aspects of 
correlations. Specifically, we present single-particle re- 
duced density matrix in the context of long range order 
in Sec.IIIA. In Sec. IIIB, We discuss the natural orbitals 
and their populations. The momentum density distribu- 
tions are shown in Sec. IIIC. In Sec. IV, we go further to 
study the two particle correlations. A summary of our 
main results is given in Sec. V. 



II. THE MODEL AND ITS EXACT SOLUTION 

We consider a system of two-component Fermi gas with 
8 interaction confined in one-dimensional box of length 
L. The model Hamiltonian reads 



JY 
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N 



* = -£ir-^5> 



dx 



h 3)i 



(1) 



l<j 



where c is the interaction strength between different spins 
of fermions and N is total number of fermions. The cou- 
pling c can be expressed in terms of the scattering length 
as c = 2/aiD, where air> is the effective one-dimensional 
scattering length. The c — and c — ► oo limits corre- 
spond to zero and infinite interactions, respectively. It 
can be tuned by the scattering length. The coupling is 
repulsive for c > and attractive for c < 0. We only 
consider the repulsive case in this work. We introduce 
another important parameter 7 that characterizes the dif- 
ferent physical regimes of one-dimensional quantum gas, 
7 = c/p where p = N/L denotes the linear density of the 
system. 

The wave function of system (1) is antisymmet- 
ric, and can be characterized by the total number 



of fermions N and the number of spin-down fermions 
M. Suppose the antisymmetric wave function is 
\l/(:ri<7i, X2CF2, ■ ■ ■ ,xn&n), where Xj and ij are the po- 
sition coordinates and spin coordinates of j'-th fermion, 
respectively. Then the static Schrodinger equation for 
the wave function ^ reads 1 2011: 



"E^ + 2c E S(x i -x j )]* = E*. (2) 

i=l l l<i<j<N 



The explicit form of the wave function of the system in 
region xq 1 < xq 2 < ... < x Qn reads [17| 

$>(xiai,x 2 (T2,-- ■ ,xn<Jn) 

J2 (- 1 ) P (- 1 ) Q ^Q 1 ,-^Q N ( r l fc Pi'--- r W fc Piv) 



P,n,...,rN 

xexpii^rjkpjX-j). 



(3) 



Here P and Q are permutations of 1, 2, . . . , N and Tj — 
±1 indicate whether the particles move right or left. The 



function A 



"■Q 1 ,...,ffQ 



(rifc Pl , . . . , r N k PN ) is given as 



^<TQ 1 ,...,£TQ JV Cl^Pl I ' ■ ' 1 r Nkp N ) 

= / J bye qi A qi , . . . , e qM A qM ) 

q,e a =±l 
M 

x Yi F r 1 k P1 ,...r N k PN (£ qa A (la ;m a ), (4) 

a=l 

with m Q being the indices of the M down spins in an 
increasing order (m a < m Q+ i), q denotes the permuta- 
tions of 1, 2, . . . , M, and ^ means the summation over 
all the permutations. The functions F are 



m-l , A / ~ H 

, . , T-r kj - A + «c/2 1 



(5) 



■y kj — A — ic/2 k m — A — ic/2 

while the coefficients S 1 are determined by the equations 
S(..., A a , A/3, . . .) _ Aj8 - A Q + ic 



S(...,A/3,A a ,...) Ap-Aa-ic 



(6) 



The hard wall boundary condition sets the requirement 
of the wave function as, 

*(0cti, • • • , x N a N ) = *(xicri, • • • , Lctjv) = 0. (7) 

The eigenvalue equation (3) can then be solved by the 
coordinate Bethe ansatz method [la]. The eigenenergy 
is given by 



N 

3=1 



(8) 



where the quasimomenta kj should satisfy the Bethe 
ansatz equations [13, Ell 
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kjL = irlj — \ arctan 
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\ arctan 
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Here j = 1,2,...,N and a = 1, 2, . . . , M, where M is 
the number of the spin-down fermions. The two sets 
of variables {kj} and {A a } are the quasi-momenta and 
spin rapidities, respectively. {Ij} and {J a } are quantum 
numbers. For the ground state, the quantum number 
configuration is described by a continuous sequence: 

Ij = j, 1 < j < N; 
J a = a, 1 < a < M. 

The values of quasi-momenta and the explicit form of 
the wave function are determined by the solutions of the 
Bethe ansatz equations (9). In the case where the in- 
teractions among the atoms are repulsive, we have two 
interesting limiting regimes. One is two-component free 
Fermi gas and the other is infinite repulsion interaction 
case which is like the case of indistinguishable spinless 
fermions. 

If c = 0, the system (1) degenerates to the noninter- 
acting two-component Fermi gas. The quasi-momentum 
takes the value of the real momentum and the fermions 
occupy the single particle momentum states according to 
the Pauli exclusion principle. The wave function of N 
noninteracting fermions with M spin-down fermions is 
given by the Slater determinant as: 



* = det 



/ M*i T) 
i>i(xi I) 
^2(0:1 1) 

Wn-m(xi T) 



1pl(xN T) 
lpl(xN I) 
4>2(X N t) 



\ 



(10) 



i/j N -m(xn T)/ 



where ip m — J j- sin(^), m — 1,2, ... ,N is the single- 
particle eigenstates of a particle in an infinite potential. 
If c — > 00, the system (1) degenerates to the spinless 
fermions. The fermions cannot feel S interaction of each 
other as the wave function must vanishes when two spin- 
less fermions touch each other, so the atoms must occupy 
the single particle momentum eigenstates with momenta 
kj = 3j- (j = 1,2, . . . , N). The anti-symmetric wave 
function for the N spinless fermions is 



* = dct 



/ipi(xi) ipi(x 2 ) 
i>3(Xl) ^3(^2) 

\lpN(Xl) 1pN{x2) 



ipi{x N )\ 

i>l{XN) 

iPn{xn)/ 



(11) 



III. SINGLE-PARTICLE DISTRIBUTION 

A. Single-particle reduced density matrix and 
off-diagonal long range order 

In quantum mechanics, the density matrix pjy = 
| , 5)(^'| plays an important role. We first consider the 
single-particle reduced density matrix 



Pi 



Tr 



23.. .N 



mm 



(12) 



where the trace means to do integrations over all the 
position coordinates and spin indices except one of them 
as we are considering the spin-half particle system. In 
coordinate representation, it is given by 

Pi(x,x') = 



N J Q ^*(x,X2, ...,xn)^{x' ,x 2 , ...,x N )dx 2 . . .dx 



N 



J \^(x 1 ,X 2 ,-,x N )\ 2 dx 1 ...dx N 

The diagonal matrix elements p(x, x) give the real 
space position density distribution of the particles, which 
are shown in Fig. [TJ Starting from noninteracting case 
c = 0, where the fermions occupy the single-particle 
momentum states, the position density is the sum of 
the N independent fermions lying in their own momen- 
tum states according to the Pauli exclusion principle. 
When interaction strength increases, the distributions 
are smoothened, and the half-width increases gradually. 
As the interaction increases further, the position density 
distribution tends to that of the spinlcss fermions. 




FIG. 1: (color online) Position density distributions of 
fermions with different interactions for N = 3, M = 1. 



quantum fluctuations. The single-particle reduced den- 
sity matrix is shown in Fig. [5] for different repulsive in- 
teractions. The single-particle reduced density matrix 
expresses the self-correlations for a single particle in real 
space. In classical mechanics, it is a Dirac delta func- 
tion 8(x — x'). So one can observe from Fig. 1 that a 
strong enhancement of p(x,x') exists along the diagonal 
line x = x' for all interaction strengths. 
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FIG. 2: Single-particle reduced density matrix p(x,x') as a 
function of interaction strength for N — 3, M — 1. Shown are 
the interactions (a) c = 0, (b) c = 5, (c) c = 10, (d) c = 50, 
(e) c = 100 and (f) c = 1000. 



B. Natural orbitals and their populations 

The spectral decomposition of the single-particle re- 
duced density matrix gives 



dx'pi(x,x')(t>i(x) = \i^i{x),i = 1,2,..., (13) 



The off-diagonal matrix elements show the correlations where 4>i{x) are the so called natural orbitals, which is 
between particles at different positions which come from obtained from the eigenfunctions of the single-particle 
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FIG. 3: (color online). Evolution of the populations of the six 
lowest natural orbitals t/>i with the interactions for the case of 
AT = 3,M= 1. 



reduced density matrix and represents an effective single- 
particle states, and Xi is population of the z-th natural 
orbital and satisfies J2i ^» = N. 

We first consider two limit cases. For noninteract- 
ing two-component Fermi gas, the single-particle wave 
functions ijji( x ) are the natural orbitals of the system 
and Xj = 2, for j = 1,2, ...,M and Xj = 1, for j = 
M + 1, ..., N — M , all the higher eigenvalues being zero 
as their reduced density matrix can be given as 



p x {x 


x') = 


M 


N-M 

)+ £ 






j=i 


j=M+l 



^(x)^(x'), (14) 



For spinless Fermi gas, the single-particle wave functions 
ipi (x) are also the natural orbitals of the system and Xj — 
1, for j = 1, 2, ..., N and all the higher eigenvalues being 
zero, 



A' 



Pi(x,x') = \ tfj*(x)tpi(x'), 



(15) 



For interacting two-component Fermi gas, we find that 
interaction leads to significant differences in the spectrum 
of natural orbitals and their eigenvalues. The popula- 
tions of lowest six natural orbitals Aj(c) with the inter- 
actions are shown in Fig. [3] We see that the first two 
natural orbitals (i = 1,2) are dominant and the other 
higher orbitals (i > 2) are separated from them. When 
the interaction increases, the populations of the lowest 
two dominant orbitals decreases while those higher or- 
bitals increase slowly. However, the lowest two orbitals 
are still dominant. We show the profiles of the two dom- 
inant natural orbitals for different repulsive interactions 



in Fig. 21 In the uncorrelated case c = 0, the system 
is in a number state, and natural orbitals coincide with 
the single-particle eigenstates. When the interaction in- 
creases, the natural orbitals will be somewhat distorted 
but the parity of the orbital is conserved. 
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FIG. 4: (color online). The profile of two dominant natural 
orbitals for different interactions for N = 3, M = 1 fermions. 



C. Momentum density distribution 

In experiments, the momentum density distribution 
detect the off-diagonal correlations of the single-particle 
reduced density matrix. The momentum density distri- 
bution is calculated from the single-particle reduced den- 
sity matrix 



n(k) = 2tt~ 



dx 



dx'p{x,x')e- lk[:L 



(16) 



From Eq. fTB")) . we see that the momentum density distri- 
bution n{k) can be understood as the Fourier transform 
of the integrated off-diagonal correlation function. Fig. [5] 
shows the evolution of the momentum density distribu- 
tion from free fermions to spinless fermions. When the 



interaction increases, the distribution becomes smoother, 
and the half-width becomes larger. When the interaction 
c extends to infinity, we find that the momentum density 
distribution is still different from the spinless fermions 
although they coincide in real space. This is to be ex- 
pected as the coordinate exchange symmetry of the two- 
component fermions is different from that of the spinless 
fermions. 

Comparing with the momentum density distribution 
of a boson gas in a hard wall trap [13], we see that 
the momentum density distributions for two-component 
fermions are smoother and the half-widths are larger. 
This is due to the Pauli exclusion principle of fermions, 
that no more than two fermions can occupy same momen- 
tum state, so in momentum space the density distribu- 
tion are broader. For infinite interaction, the hard-core 
bosons have the same position density as the spinless 
fermions, but the momentum density distributions are 
different [13, 1221- The reason is that their wave func- 
tions have different symmetry, that is, the wave function 
of the bosons is symmetric while that of the fermions is 
antisymmetric. 
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FIG. 5: (color online). Momentum distribution of fermions 
with different repulsive interactions c for N — 3,M — 1. 



IV. TWO-PARTICLE DENSITY DISTRIBUTION 

Now, we consider the two-particle correlations. 
The two-particle reduced density matrix is p 2 = 
£?*34...iv-|\E r )(^ r |- In coordinate representation, the two- 
particle density, which is the diagonal of the two-particle 
reduced density matrix, also called the pair distributions, 
is obtained from the many-body wave function ^ as: 

p 2 {x u x 2 ) = 

N J Q <S>*(xi,x 2 , ...,Xn)^(xi,x 2 , ...,x N )dx 3 . . .dx N 



The two-particle density expresses the joint probability 
of finding one particle at position x\ and any second one 
at x 2 . We show the two-particle density in Fig. [6] for dif- 
ferent interactions. We find that the two-particle density 
is minimized along the diagonal line for all the interac- 
tions due to the Pauli exclusion principle. When the in- 
teraction increases, the distributions split into two parts, 
which demonstrates that the particles are more localized 
in the well. 
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FIG. 6: Two-body density p(x\,X2) as a function of interac- 
tion strength for three fermions with one down-spin. Shown 
are the interactions (a) c = 0, (b) c = 1, (c) c = 10, (d) 
c = 50, (e) c = 100, and (f) c = 1000. 



/ \^(x 1 ,x 2 ,...,x N )\ 2 dxi 
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V. CONCLUSION 

We have investigated the ground state properties of the 
one-dimensional two-component Fermi gas in an infinite 
potential well. From the solutions of the Bethe ansatz 
equations, we obtained the explicit form of the wave func- 
tion of the system for a few atoms. Then we studied 
the single-particle reduced density matrix, momentum 
density distribution and two-particle density matrix for 



different interactions. We aiso found that the momen- 
tum density distributions of the strongly interacting two- 
component fermions were very different from that of the 
spinlcss fermions, while the position density distributions 
of both were the same. 
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